Electron-vibration energy exchange models in nitrogen-containing plasma flows 
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The physics of vibrational kinetics in nitrogen-containing plasma produced by colli- 
sions with electrons is studied on the basis of recently derived cross sections and rate 
coefficients for the resonant vibrational-excitation by electron- impact. The tempo- 



derivation of reduced models and to the accuracy of a relaxation time formalism. A 
analytical fit of the vibrational relaxation time is given. 



^ ral relaxation of the vibrational energy and of the vibrational distribution function 

is analyzed in a state-to-state approach. The electron and vibrational temperature 

are varied in the range of [0-50000] K. Conclusions are drawn with respect to the 
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I. INTRODUCTION 

Modeling of high speed flows of reactive gas mixtures is required for the description of 
space vehicles entries in planetary atmospheres, as well as for the interpretation of ground 
based experiments in high enthalpy wind tunnels and shock tubes. 

It has long been recognized that vibrational energy relaxation, beside its obvious role in 
the global energy balance, plays a major role in the chemical kinetics of non-equilibrium 
zonesi. Among the many processes that influence the molecular vibrational energy, in this 
paper we focus on inelastic collisions with free electrons, the so-called electron-vibration 
(e-V) processes. Even in conditions of weak ionization, as they occur in Earth entries at 
orbital speed or in Mars entries, excitation by electron-impact is very efficient a mechanism 
for populating vibrationally excited molecules. In addition, the aerospace community shows 
a growing interest for flows with sensibly higher ionization degrees. This is due to the 
interest in superorbital speed Earth entries as they happen in the return phase from Lunar 
and Martian exploration missions^-, and to the investigation of plasma flow control concepts 
(see Ref. |4j for a review). 

For nitrogen-containing plasma flows, the direct process of N2 vibrational-excitation by 
electron scattering is a optically forbidden transition, due to the lack of dipole momentum 
for the N 2 molecule. For incident electron energies lower than 4 eV, however, the vibrational 
transition occurs via a resonant process involving the 2 U g state of the N^ ion: 

e~ + N 2 (X X E+; v) — ► N 2 "( 2 n 9 ) — ► e" + N 2 (X *£+; w) , (1) 

where v and w are the initial and final vibrational level of N 2 . By means of the resonant 
process in ([1]) the vibrational-excitation cross section is enhanced by several orders of mag- 
nitude. This process is very efficient and makes the vibrational temperature equilibrate 
very rapidly with the free electron temperature^ to the extent that many simple models 
assume the N 2 vibrational temperature instantaneously equilibrates with the free electron 
temperature. A first serious theoretical study for e-V relaxation by electron-impact has been 
proposed by Lee^. More refined models have been developed by Mertens^ and by Bourdon 
and Vervisch^. 

The experimental data on resonant vibrational-excitation cross sections and rate coeffi- 
cients are difficult to obtain and, in the literature, they are limited to the first few levels^. 



In order to build a high temperature model information on the upper vibrational levels 
is mandatory. Recently, theoretical calculations have been performed for electron-nitrogen 
resonant scattering considering all vibrational transitions up to the dissociation limit in the 
energy range — 10 eV— that show very good agreement with available experimental data. 
The new theoretical results are used in this paper to set up a kinetic study of electron-impact 
vibrational energy relaxation in nitrogen-containing plasma with the purpose of assessing 
the accuracy of existing models. 

The paper is organized as follows: in section [III the basic equations of the kinetic model 
are presented; in section II III the full set of rate constants for all vibrational transitions are 
introduced and used to assess modeling assumptions commonly used in the literature. The 
solution of the full kinetic model is analyzed in section IIVI together with a reduced multi- 
quantum vibrational transition approach. The relaxation time formalism is introduced in 
section |V]and the results of its application to the vibrational energy relaxation are compared 
to those obtained within the full kinetic model. Section IVT1 summarizes the main conclusions 
and perspectives. 

II. MODELING OF ELECTRON- VIBRATION RELAXATION 

In order to model the e-V kinetics a nitrogen gas in contact with a free electron bath 
is considered. The electrons are considered at equilibrium at the constant temperature T e 
and the N2 molecules are supposed to be in the electronic ground state X 1 EJ", with a initial 
vibrational distribution at temperature T v . Due to the very efficient resonant mechanism, 
the electron-N 2 collisions drive the vibrational distribution towards the equilibrium at the 
electron temperature T e . Since the electron affinity of the N atom is negative, the resonant 
dissociative electron attachment for electron-^ is suppressed^. Therefore, the effects of 
molecular dissociation on vibrational relaxation are not taken into account. 

The time evolution of the molecular vibrational distribution function (VDF) is described 
by the following equations: 



n e > [ JXyuy n w j\ vw n v ^ , yz j 



dt 

w 



where n e is the electron density and K V;W (T e ) is the rate constant for the v —¥ w resonant 



vibrational-excitation (RVE) process: 

K v , w (T e ) = ^=(k B T e )~ 3/2 J ea v , w (e)e-^ T ^de, (3) 

where a v ^ w is the corresponding cross section as a function of the electron energy e. In 
Eq. 02]) the sum extends over all 67 vibrational levels of N2. The detailed balance principle 
imposes the following condition among the direct and inverse processes: 

where n* v is the equilibrium population of the vibrational level v at temperature T e : 

n* v ~ e" e " /AsTe , (5) 

e v being the energy of the N 2 vibrational level v and k^ the Boltzmann constant. The total 
N2 vibrational energy is readily obtained from the VDF: 

Kib = ^2,n v e v . (6) 

V 

III. RATE COEFFICIENTS 

In this section the rate coefficients K VtW (T e ) that enter the kinetic equations (Cp and ([6]) 
for the resonant process in Eq. 01} are presented. They are taken from Ref. |_10| (paper I 
in the following) where full details on the calculations and the validation of the model by 
comparison with other approaches can be found. The method used in I to calculate the 
electron-N2 resonant scattering is the well-known 'boomerang model', an approximation 
of the nonlocal-complex-potential modet^. Such a method has been used to study many 
different systems and a complete description together with the application limits can be 
found in the papers^— and references therein. Molecular rotations have been taken into 
account by adding the centrifugal barrier to the potential, parameterized by the rotational 
quantum number J and transitions with A J = have been considered. Results show that 
the rate constants for vibrational-excitation do not depend strongly on rotational excitation 
so, in the following, only J = is considered. 

In I, a full set of RVE cross sections among all vibrational levels of N 2 has been calculated. 
The corresponding rate coefficients have been obtained by averaging each cross section over 
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FIG. 1. Theoretical rate coefficients for electron-N2 resonant vibrational excitations as a function of 
electron temperature as obtained in Ref. lid . Left: Mono-quantum transitions; Right: n-quantum 



transitions starting from the ground vibrational level. 



a equilibrium electron energy distribution as defined in Eq. ([3]) in the range [0-50000] K of 
electron temperature. Figure [1] shows typical rate constants for mono-quantum transitions 
(v — > v + 1) and for multi-quantum transitions starting from the ground vibrational level 
(0^0 + n). 

Since a large number of transitions are involved in the full data set, it is natural to look 
for regularities that allow using scaling relations. On the other hand, in the past both 
theoretical and experimental studies on RVE by electron-impact were, as a rule, limited to a 
small subset of all possible processes: that required the introduction of scaling-laws among 
the rate coefficients in order to obtain the missing data. Mertens^ proposed to assume 
a linear relation between \og(K 0iV ) and the final vibrational level v. Bourdon et a/.- have 
shown, using the Allan datar^, that this rule is approximately valid up to v = 7 at T e = 1 eV. 
In Fig. [2] we report the calculated rate coefficients K 0tV as a function of the final vibrational 
level v for different electron temperatures T e , together with a linear extrapolation obtained 
using the first 5 points. A linear relation holds with good accuracy for levels v ^ 5. At higher 
electron temperatures this limit can be raised to v ~ 10 but transitions with larger jumps 
are overestimated by many orders of magnitude. As we shall see in Sec. IIV| transitions with 
quantum jumps up to Av max = 20 have to be considered for a accurate estimation of the 
vibrational energy relaxation. 



For transitions starting from vibrational levels other than the ground state, GordietsA^ 
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FIG. 2. Rate coefficients Kq v as a function of the final vibrational level v calculated at different 
electron temperatures. Dashed lines are linear extrapolations obtained using the first 5 points. 




FIG. 3. Ko 7 n/K VtV + n as a function of the vibrational level v, for different values of the electron 
temperature T e , for mono-quantum (left) and bi-quantum (right) transitions. 



proposed the following scaling law: 



K, 



*>-v,v+n 1 , ' \ ' / 

1 + av 
a being an adjustable parameter. In Fig. [3] the ratio K^ n / K v>v+n as a function of the initial 
vibrational level v, for different values of the electron temperature and for different values of 
the vibrational jump n, is reported. The results show that the linear relation, Eq. ([7|), holds 
approximately only in limited regions of the parameters' space and that the a parameter is 
temperature-dependent . 

The scaling-laws proposed by Gordiets and Mertens are heuristic guesses coming from 
gas discharge physics^ where typical conditions are electron temperature around 1 eV and 
the gas temperature below 1000 K. Figures [2] and [3] show that indeed linear scaling-laws 
are reasonable assumptions when the electron temperature is in the range where the rate 
coefficients have a broad maximum (10000 < T e < 30000 K see Fig. [TJ and for vibrational 
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levels up to v = 10. For high speed flows applications, however, models should be valid in 
a much broader range of the parameters' space. The linearity deteriorates for low electron 
temperatures and high vibrational levels and extrapolation of these relations can lead to 
severe overestimation of the rate constants. In the following, only the rate coefficients from 
paper I are used. 

IV. VIBRATIONAL ENERGY RELAXATION 

Using the rate coefficients presented in section IHff the time evolution of the VDF of N 2 
molecules and of the vibrational energy E v n, is studied by solving the coupled rate equations 
in Eqs. (0). In particular the role of the multi-quantum transitions is investigated. The 
initial condition for the VDF is supposed to be a Boltzmann distribution at temperature 
T v . At given electron temperature, the rate constants for transitions involving the exchange 
of many vibrational quanta are much smaller than the others due to the large energy gap: 
see for example Fig. [TJ It is therefore natural to ask wether it is possible to set a upper 
limit to the number of quanta exchanged in order to reduce the amount of data required for 
the modeling. Figures H] and [5] show the 'full' solution, i. e. obtained including the full set 
of multi-quantum transitions, and the effect of considering a restricted range of vibrational 
transitions in the vibrational kinetics on the VDF and on the corresponding E v n, for two 
typical conditions of heating (T v = 5000 K and T e = 30000 K) and cooling (T v = 20000 K 
and T e = 2000 K). VDF is normalized to the population of the vibrational ground state and 
vibrational energy E vi b is normalized to its equilibrium value E* ib . Times are normalized 
by: 

roi(T e )= l , (8) 

so that the plots are independent of electron and molecular number densities. Figure H] 
shows the relaxation of the VDF. The latter, although a Boltzmann distribution at the 
beginning (by construction) and at the end (by the equilibrium condition) of the relaxation, 
shows non-equilibrium character both in heating and cooling conditions. In both cases, high 
vibrational levels relax more slowly than low ones and the VDF exhibits strong depletion 
(overpopulation) of high levels with respect to a Boltzmann distribution at the same total 
energy. Large multi-quantum jumps have a strong influence on this relaxation and restricting 
the maximum number of exchanged vibrational quanta overestimates the degree of non- 
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FIG. 4. Temporal evolution of normalized vibrational distributions function as obtained by different 
assumptions on the maximum number of vibrational quanta exchanged in collisions. Full line: 
all transitions allowed; dashed line: Av ma x = 10; dot-dashed line: Av max = 20. Left panel: 
T e = 30000 K; T v = 5000 K. Right panel: T e = 2000 K; T v = 20000 K. 
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FIG. 5. Temporal evolution of the normalized vibrational energy obtained by different assumptions 
on the maximum number of vibrational quanta exchanged in collisions. Full line: all transitions 
included; dashed line: Av max = 10 levels; dot-dashed line: Av max = 20 levels. Left panel: 
T e = 30000 K; T v = 5000 K. Right panel: T e = 2000 K; T v = 20000 K. 

equilibrium by many orders of magnitude. The plot on the right in Fig. H] shows that this 
error can involve vibrational levels as low as v ~ 10. It is interesting to note that a fairly 
good agreement is found in Fig. [5] for vibrational energy, even if significant discrepancies are 
observed on the VDFs, as vibrational energy is mostly due to low vibrational levels. 

Figure M reports the maximum error, calculated over all the temporal evolution of E vi b, 
between the results obtained by considering the full set of rate coefficients and those obtained 
including only transitions with Av max = 10 and Av max = 20 in the kinetic model. Results 
show Av max = 10 is a reasonable approximation with a maximum error smaller than 10% in 
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FIG. 6. Maximum error during the temporal evolution of E v n, obtained by considering multi-quanta 
transitions with Av max = 10 (left) and Av max = 20 (right). 



a wide region of the parameters' space. Large errors appear for cooling conditions (T v > T e ) 
when the starting vibrational temperature is large enough. Although so large vibrational 
temperatures are difficult to realize in practice, this is an indication that very large multi- 
quantum jumps are involved in the relaxation of strongly populated high vibrational levels. 
With Av max = 20 the maximum error stays below 1% for most conditions of interest. 

The vibrational energy relaxation is therefore described with reasonable accuracy by 
taking into account only transitions where less than 10 vibrational quanta are exchanged; if a 
accurate description of the vibrational distribution is required, however, at least Av max = 20 
must be considered. 

To conclude this section, Fig. [7] summarizes the temporal evolution of the vibrational 
energy, as obtained by the full kinetic model, for two electron temperatures 3000 K and 
12000 K and for initial vibrational temperature in the range T v = to 50000 K to model 
different heating and cooling conditions. It is clear that the rate of relaxation is mainly 
determined by the electron temperature. But the initial vibrational distribution also plays 
a role since different initial conditions relax along different paths. 
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FIG. 7. Time evolution of the normalized vibrational energy, E v ij,/E* ib , for two different values of 
the electron temperature T e , parameterized on the N2 initial vibrational temperature T v . 

V. RELAXATION TIME FORMALISM 

The standard way of modeling the vibrational energy relaxation is to describe it by a 
Landau- Teller (LT) rate equation: 



dE„ 



vib 



E 



vib 



E 



vib 



(It 



(9) 



where E* ib is the equilibrium vibrational energy at the electron temperature T e , and r e is 
the e-V relaxation time. Early theoretical work of Lee^ studied the relaxation from a T v = 
condition assuming a harmonic oscillator model for the N 2 molecule and using the cross 



sections from Ref. 



for r e : 
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Results can be summarized with the following analytical expression 



Pe^ 



Lee 



^U(|) 2 E J ^(T e ). 2 | 



-1 



(10) 



where p e = nekgTe is the electron partial pressure and 8 V is the characteristic vibrational 
temperature of the N 2 molecule. The relaxation time defined in Eq. ( ITU]) depends only on T e 
and it does not take into account vibrational excitation. Lee's formulation has been tested 
and improved in the works of Mertens^ and Bourdon et alM^-. In these works, cross sections 
have been taken from the theoretical work by Huo et a/.— or from the experimental data 
of Allan^. The Huo et al. results are ab-initio calculations for — )■ v ^ 5 transitions 
with rotational quantum numbers J = 0, 50, 150. The corresponding rate coefficients are 
calculated in the electron temperature range 0.1 — 5.0 eV. In addition, results have been 
obtained for J = 50 and initial vibrational quantum number up to v = 12 with changes in 
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FIG. 8. Vibrational relaxation time as a function of electron temperature for different values of 
the initial vibrational temperature T v . The values are obtained by solving the full set of kinetics 
equations. 



vibrational quantum number of —5 and +5. In particular, the results obtained for J = 50 
have been used in the modeling calculations^. The Allan's results are experimental values 
for —?• v transitions up to v = 13 obtained at room temperature. In both cases, cross 
sections for all other transitions must be obtained by assuming a scaling-law^ 1 ^. As a result, 
analytical expressions have been derived for the relaxation time in Eq. (Q that depends both 
on T e and T^~-. 

Analyzing the full set of RVE rate coefficients, we concluded in Sec. IHII that such scaling- 
laws have strong limitations. It is therefore interesting to investigate the limits of a LT 
description of the vibrational energy relaxation and to provide, where possible, updated 
values for the relaxation time. To this end, time histories of vibrational energy for differ- 
ent conditions of electron temperature T e and initial vibrational temperature T v have been 
calculated with the full kinetic model. As it will shown in the following, in general, the 
relaxation is not a linear process. Therefore some arbitrariness is implicit in the definition 
of the relaxation time. In this paper the value of r e was obtained in each case by linearizing 
the relaxation rate around its starting point, in order to minimize the discrepancies between 
the LT and the full model values for the vibrational energy. 

The final results are reported in Fig. [8] and in particular in Fig. [9] they are compared with 
those of Lee^ and Bourdon et a/.- for T v = 0. Present results lie between the other two, 
somewhat closer to the more recent results of Bourdon et al. 

Due to the resonant character of the cross sections, the relaxation time shows a mini- 
mum. In Table [J the point of the minimum (T e mm , p e r™ in ) of the relaxation time, for some 

11 





14 


1. 


: 




12 


! 


: 










7 


10 


1 ! 




5 




1 i 


,..• ■ 










"i 1 


" 


I » 




~ 


















h 


\ ' 


S 














\ \ 


*"** ^^ 




4 




-"'"'' ^-^^^ 




2 













5000 10000 15000 

Electron temperature (K) 



FIG. 9. Vibrational relaxation time as a function of the electron temperature with initial vibrational 

P 

temperature T v = K. Full line: this work; dashed line: results from Ref. |6fc dotted line: results 
from Ref. [9|. 
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FIG. 10. Best fit of T™ m and p e r™ m as a function of the initial vibrational temperature, as given 
in Table ffl 



vibrational temperatures, is shown. These points can be fitted by the following laws: 

rrrniinfrri •• 



T v + a 

Pe rf n (T v )=c + dT v 



where the coefficients take the following values: 



(11) 
(12) 



a = 22159.90 K, 
b= 1.60 x 10 8 K 2 , 
c= 1.20 x 10" 5 Jm~ 3 s, 
d = 5.9 x 10~ 10 Jm _3 sK _1 . 

Figure [TO] shows the best fit superposed to the calculated points. 



(13) 
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T v (K) T e min (K) p e Tf n (1(T 5 Jm- 3 s) 






7055.64 


1.29 


1000 


7010.24 


1.29 


2000 


6776.13 


1.29 


5000 


5968.45 


1.35 


10000 


4932.81 


1.65 


20000 


3800.00 


2.62 


50000 


2189.37 


4.11 



TABLE I. Minimum of relaxation time (T™ in , p e r™ m ), for some initial vibrational temperatures 



T v . 



For electron temperature T e ^ j 1 ™ 111 the relaxation time curves in Fig. [8] can be fitted, 
with 1% accuracy, by the following polynomial expression: 

3 
p e r e (Te, T v ) = Pe rr(T v ) + J2 C *( T -) i T e - Tf n (T v )Y , T e ^ T e min (14) 
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(15) 

The LT equation, Eq. Qj, can be integrated analytically to give the vibrational energy 
relaxation. For each set of initial conditions (T e ,T v ), a value for the relaxation time is 
calculated and kept constant during all relaxation; the resulting vibrational energy and the 
corresponding value obtained by solving the full set of kinetic equations were compared 
and the maximum deviation was estimated. Figure [11] show the results of this analysis 
as a contour map. A single relaxation time, as in LT approach, is able to describe the 
e-V process both in relaxing and exciting conditions when the non-equilibrium is not very 
strong (T e ps T v ). The range of validity increases for larger T e . Under strong non-equilibrium 
conditions, however, large discrepancies (larger than 10%) are seen. 

In order to explain the results of Fig. [Til the relaxation process was analyzed in greater 
detail for three cases. Table [TTJ summarizes the conditions for the cases studied, the value for 
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FIG. 11. Maximum deviation on the vibrational energy as obtained with a Landau- Teller approach 
compared with that obtained by solving the full vibrational kinetic model. 



the relaxation time and the computed error with respect to the exact solution. Figure 
shows the time evolution of the normalized vibrational energy for two heating conditions 
with T v = K and T e = 5000 and 15000 K and a cooling condition with T v = 8000 K and 
T e = 2000 K. We note that the relaxation predicted by the LT equation is faster than the 
exact solution both in exciting and relaxing conditions. Figure [13] shows the corresponding 
normalized VDF at different times during the relaxation. For the LT model, Boltzmann 
distributions at all times during the relaxation are assumed with vibrational temperature 
calculated from vibrational energy. In exciting conditions, at the beginning, high vibrational 
levels (v > 1) relax faster than low levels and get strongly overpopulated. In this conditions, 
the LT predictions agree well with the 'full' result since the low- lying levels make the largest 
contribution to E v n,. At later stages, the reverse is true: again the LT results agree with the 
'full' results for the first two levels, but now all the rest of the VDF is overestimated. This 
is what makes the LT value for E v n larger. Note that also for case 1, where the relaxation 
predicted by the LT equation agrees well (within 2%) with the 'full' result, the VDF are 
strongly non-equilibrium during most of the relaxation. This can have a strong influence 
when a vibrationally favored chemical reaction (e.g. dissociation) is considered. In relaxing 
conditions, instead, the high-energy levels relax slower than low levels. As a result, the 
LT relaxation is faster than the 'full' value. Again, strongly non-equilibrium VDF appear 
during most of the relaxation. 

14 



case 


T e (K) 


T v (K) 


Te (s) 


^max 


1 


5000 





0.023 


0.01 


2 


15000 





0.013 


0.26 


3 


2000 


8000 


0.14 


0.62 



TABLE II. Definition of the cases treated in the text. r e is calculated by linearizing the energy 
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FIG. 12. Time relaxation of normalized vibrational energy for the cases described in Table HH Full 
line: exact solution as obtained by solving the full set of kinetic equations; dashed line: constant 
relaxation time in a Landau- Teller approach. 



This is the origin of the failure of a relaxation time description. Obviously, the modeling 
done in Refs. |6| has been improved in Refs. |7|j9| by deriving a relaxation time that depends 
both on T e and T v . 

Based on the results presented in this work for the relaxation of the VDF, however, we 
do not expect the inclusions of the dependence of the relaxation time on T v to significantly 
improve the predicting capabilities of the LT model. Indeed, in Fig. [7] we show the E v n, 
relaxation for two values of T e and for several values of T v . It is apparent that these quantities 
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FIG. 13. Relaxation of the N2 vibrational distribution function, normalized to the ground vibra- 
tional state, for the cases described in Table HH Full line: exact solution as obtained by solving the 
full set of kinetic equations; dashed line: constant relaxation time in a Landau- Teller approach. 

alone are not able to describe the system relaxation towards equilibrium. 

VI. CONCLUDING REMARKS 



Recent calculations of RVE cross sections in nitrogen covering all vibrational transitions^ 
up to the dissociation limit have been used to study the vibrational energy relaxation as 
produced by these resonant processes. Analysis of the full set of rate constants shows that 
approximate scaling-laws adopted in the past have a limited range of validity. 

Also, the possibility to model the relaxation by a linear rate equation with a single 
relaxation time has been investigated. The results from a relaxation time formalism and 
from the full kinetic model are compared to conclude that the former may be used with 
reasonable accuracy (~10%) to model the vibrational energy relaxation in most cases of 
practical interest. A analytical fit of the calculated relaxation time, that depends both on 
free electron and vibrational temperatures, has been provided in Eqs. (!TT1) - (TT5|) that could 
be of interest for the modeling of this important relaxation process. 

State-to-state kinetic modeling, however, predicts strongly non-equilibrium vibrational 
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distributions during the relaxation both in heating and cooling conditions. Although not 
considered here, this can have a strong influence on the chemical kinetics of reactive plasma 
flows, in particular where vibrationally favored reactions (e.g. dissociation) play a role. 

The next step will be the inclusion, in the relaxation kinetic of molecular rotation and of 
the dissociation reaction. 

Finally, it is worth mentioning that the original cross section data allow the modeling of 
RVE processes also in conditions where the electron energy distribution function deviates 
significantly from equilibrium. For the present study only rate coefficients have been used. 
These are obtained from the cross sections by averaging with a Maxwell distribution for the 
electron velocities. 
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